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Abstract. We present a detailed model for the formation 
of massive objects at the centers of galaxies. The effects 
of supernovae heating and the conditions of gas loss are 
revisited. The escape time of the gas is compared with 
the cooling time, which provides an additional condition 
not previously considered. Its consequences for the allowed 
mass range of the halo is calculated and parameterized in 
terms of the spin parameter, A^,, the redshift of collapse, 
Zc, the fraction of baryons in stars, /*, and the efficiency 
of supernovae, v. It is shown that sufficient gas is retained 
to form massive dark objects and quasars even for moder- 
ately massive halos but a decline is expected at low red- 
i; shifts. Subsequently, a gaseous disk forms with a radial 
^ ■ extent of a kpc, spun up by tidal torques and magnetized 
O , by supernovae fields with fields strengths of 10—100// G. In 
' a model of a self-similar accretion flow in an initially dom- 
QQ inant halo, it is shown that for typical halo parameters, 
p\| about IO^Mq accretes via small magnetic stresses (or al- 
, ternatively by self-gravity induced instability or by alpha 
viscosity) in 10* years into a compact region. A model of 
a self-gravitating evolution of a compact magnetized disk 
Qv^ ' ('"oJ; 100 pc), which is relevant when a significant fraction 
, of the disk mass falls in, is presented, and it has a rapid col- 
' lapse time scale of a million years. The two disk solutions, 
one for accretion in an imposed halo potential and the 
other for self-gravitating disk, obtained here, have general 
utility and can be adapted to other contexts like proto- 
■ stellar disks as well. Implications of this work for dwarf 
galaxy formation, and a residual large scale seed field are 
O also breifly discussed. 

^ Key words: Accretion, Magnetic Fields, Galax- 
ies:Formation, Cosmology: Theory 
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1. Introduction 



There seems to be increasing evidence that supermas- 
sive black holes arc at the centers of galaxies. Dynami- 
cal searches indicate the existence of massive dark objects 
(MDOs) in eight systems and their masses range from 
IqS _ 10^-5 Mq (Kormcndy & Richstonc 1995). Although 
this study does not confirm that the central objects are 
supermassive black holes, it has been inferred that the cen- 
tral mass is contained within 10^ Schwarzchild radii. On 
an average, the black hole mass is a fraction, 10~^ — 10~^, 
of the total mass of the galaxy and of order 10~^'^ of the 
bulge mass (Wandcl 1999). Recent observations show a 
strong correlation between the black hole mass, ni/j, from 
stellar dynamical estimates, and the velocity dispersion 
of the host bulges (in;i oc a"; where a is reported to be 
in the range 3.5-5; eg. Ferrarese & Merritt, 2000). This 
has been supported by reverberation mapping studies of 
the broad line region (Gcbhardt et al. 2000). The black 
hole masses in active galaxies as inferred from their lu- 
minosities, assuming reasonable efficiencies are also in the 
same range, 10^~^'^Mq. Arguments based on time vari- 
ability, relativistic jets and other circumstantial evidence 
indicate that they are relativistically compact (Blandford 
& Rees 1992). Specific examples include a ^ 20 pc disk 
spining at 500 km s~^ which implies a 2 x IO^Mq black 
hole in M87 (Ford et al, 1994) and evidence of IO'^Mq 
mass in a region of 0.1 pc extent in the case of NGC 4258 
(Miyoshi et al. 1995). One remarkable fact is that there is 
a decline in the quasar population between z = 2 and the 
present epoch. The presence of quasars at high redshifts 
tells us that galaxy formation had proceeded far enough 
for supermassive black holes to form in the standard pic- 
ture (Rees 1984). A detailed model of formation of these 
objects, such as the one attempted here, should address 
the issues of supernovae feedback from star formation and 
the mechanism of efficient angular momentum transport 
in order to explain the massive active nuclei as early as 
^; = 5. In the case of MDOs, there is a need to explain 
the compact sizes of 10 — 100 pc that are implied from 
dynamical studies. 



Broadly, the two main routes to the formation of 
the massive central objects that have been proposed are 
through instabilities in a relativistic stellar cluster or gas 
dynamical schemes which may involve a direct collapse of 
a primordial gas cloud or accretion of a collapsed gaseous 
disk. The main drawback of the stellar cluster models is 
that one must assume the existence of dense and mas- 
sive cluster at the outset; the angular momentum trans- 
port problem to arrive at this initial scenario is difficult to 
overcome. One gas dynamical scheme proposed by Shlos- 
man, Begelman & Frank (1990) involves accretion of gas 
through stellar bars which have been induced either by 
self-gravity or galaxy interactions) driving the gas from 10 
kpc to about a kpc size 'disk of clouds' which further ac- 
cretes by viscous dissipation due to cloud-cloud collisions. 
N-body simulations (Sellwood & Moore 1998) have shown 
that the bar weakens substantially after a few percent of 
the disk mass accumulates in the center. This mechanism 
may not last long enough to drive suffident matter into 
a compact region as the bar instability is suppressed by 
the bulge at the inner Linblad resonance. Disk accretion 
due to self-gravity and magnetic fields may then be better 
candidates to transport the mass to about a 100 pc size 
compact region. 

Loeb & Rasio (1994) considered the possibility that 

massive black holes form directly during the intial col- 
lapse of the protogalaxies at high redshifts and per- 
formed smoothed particle hydro dynamical simulations of 
gas clouds. They find that inital collapse of a protogalac- 
tic cloud leads to the formation of a rotatinally supported 
thin disk. They argue that if the viscous transport time in 
the disk is small compared to the cooling time, then the 
disk could collapse to a supermassive star, else it could 
form a supermassive disk; they do not provide gas dy- 
namical collapse model after the formation of a disk or 
star. They found that the gas fragments into small dense 
clumps (presumably are converted into stars). The impor- 
tant consequences of supernovae feedback is not consid- 
ered. Eisentein & Loeb (1995) suggest that quasars may 
be associated with rare systems that acquire low values of 
angular momentum from tidal torques during the cosmo- 
logical collapse. Typically these are IO^Mq objects and 
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form at redshifts z'> 10. Since the viscous evolution times 
in the coUapsed disks of these objects are comparable to 
star formation times, it is important to consider the effects 
of supernovae feedback in detail that could disrupt these 
objects. 

Natarajan (1999), considers the effect of feedback and 
argues that the fraction of the gas retained is proportional 
to ratio of the supernovae heat input to the binding energy 
of the gas in the halo. This assumption has interesting con- 
sequences for the Tully-Fisher relationship. The criteria 
for mass loss was considered by (Mac-Low & Ferrara 1999; 
Ferrara & Tolstoy 2000) in the context of a blow-out or 
a blow-away (when the mass is completely expelled) from 
disk galaxy with an isothermal atmosphere. The condition 
for the gas to be retained is that the blow-out velocity be 
less than the escape velocity. However, the cooling of the 
hot gas is not considered. In this paper, we consider the 
possibility that gas cools before it escapes the halo. 

Silk & Rees (1998) propose that ~ IO^Mq clouds at 
high redshifts undergo direct collapse in a hierarchical 
(bottom-up) CDM cosmology without undergoing frag- 
mentation and star formation arguing that the conditions 
in primordial clouds differ from conventional molecular 
clouds. The angular momentum in this picture is shed 
by non-axisymmetric gravitational instabilities. The more 
massive black holes ( IO^Mq) ejects mass from the host 
galaxies by assumed spherical quasar winds. In the lim- 
iting case, it provides a relationship between the mass of 
the central hole and the host galaxy mass. The time be- 
tween halo virialization and the birth of quasars is short 
compared to the cosmological timescale. Haehnelt & Rees 
(1993) proposed a scenario in which a disk forms (after 
turnaround and collapse) and loses angular momentum 
due to OL viscosity in an estimated timescale of 10^ yrs. 

Here, we discuss a detailed physical model for the for- 
mation of protoquasars (or MDOs) from a magnetized ac- 
cretion of a collapsed disk, the properties of which are ob- 
tained taking into account supernovae feedback in a viri- 
alized halo. There is observational evidence that consider- 
able fragmentation precedes quasar activity and the broad 
emission lines in quasar spectra indicate high metallicity 
(Hamann & Ferland 1992). We assume, therefore, signif- 
icant star formation and supernovae activity occurs after 
the cloud, which is spun up by tidal torques, contracts 
to a radius where self-gravity is significant. The paper is 
composed of the following parts (See Fig |l|)- 

1. The formation of a gaseous disk with a radial extent of 
about a kpc, in a host galaxy as limited by supernovae 
feed back. We investigate in the range in halo mass 
for a given redshift that still retains the hot gas. The 
effect of the evolution of gas on the coUisionless dark 
matter system is neglected. 

2. In prev ious work, gravitational instabilities (examined 
in §0) in the disk was considered as the main source 
of viscosity. In § 4.1 , justification is made for a magnetic 
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Fig. 1. The evolution of a 10^°Mq cloud leading to a for- 
mation of a protoquasar. In first step, this cloud cools 
faster than the dynamical time, tc < ta- After star forma- 
tion, one of the conditions is satisfied- that the supernovae 
heat energy is less than the binding energy, Esn < Ebe 
or that the cooling time of the hot gas is less than the 
time to escape the halo, tg <te- 



viscosity from supernovae fields and the estimated ac- 
cretion rate turns out to be significant; also the large 
scale field strength derived from the effective seed of 
small scale fields is used to explain observations (§^. 
The collapse of the disk is calculated with a gener- 
alized viscosity prescription (which includes the indi- 
vidual cases of magnetic, a and self-gravity induced 
instabilities, §^) under a halo dominated gravitational 
potential (§^ into a compact central region at rapid 
rate of about a M© yr~^. A self-gravitating magne- 
tized disk solution for this central object that collapses 
in 10® yrs, is presented in §^ There is summary and 
discussion in §0 and conclusions in 
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2. The virialized spherical halo and formation of 
the gaseous disk 

We assume a standard spherical model for the formation 
of a virialized spherical halo to begin with, in which the 
gas cools to form the disk. Subsequently, we include in 
the calculations the supernovae heating and consider the 
conditions under which sufficient gas is retained to form 
the disk. 

A particularly simple and useful version of the spher- 
ical model below assumes that the matter distribution 
is symmetric about a point and is a pressureless fluid. 
The shell enclosing mass of the overdense region, Ai, ini- 
tially expands with the background universe, slows down, 
reaches a maximum radius before it turns around, and col- 
lapses. The collapse proceeds until a time when it reaches 
virial equilibrium. For an average density contrast, 6, and 
using the fact that the background density, pb cx for 
a flat cosmological model, one can make the following 
estimates of the typical parameters of the collapsed ob- 
ject (Padmanabhan & Subramanian 1992, Padmanabhan 
1993, Peebles 1980) 
po = 9.1 X 10-^° n h^Q g ciR-^ 

To = 0.94 X lO^^/iTO yr 

r„ = rt/2 - {78/do) X}^' h-^^\pc (1) 
Te = 2.19 ro/(5o/' 

So = (37r/2)3/2 ^ ^ ^ ggg ^ ^ -,^_Qg2 (1 -h zt), 
where po and tq are the current values of the comoving 
density and time, the subscripts t and c indicate turn 
around and collapse values, and is the radius of virial- 
ization. For a given model of the cosmological evolution 
of the initial density perturbation, one obtains values for 
(5o or alternatively one can specify the collapse redshift. 
As an example, zt = 6.3 corresponds to Zc = 3.6 and a 
(5o =7.8. We take Jl;, = 0.1 K^^ and the mass in baryons, 
Mg, in units of lO^M©, is related to the total mass by 

Ms = M ^ A M^o. 
The baryonic mass, A/g, will fuel a black hole, after it has 
formed, at a rate limited by the Eddington luminosity. If 
the accretion proceeds at a tenth of this rate, then the 
luminosity in units of lO^^ergss"^ is 

L45 = 1.3 m;,8, (3) 

where Ill^ts is the mass of the black hole in units of 10^ Mq . 



2.1. Calculation of the collapse factor ry / 

The gas in the massive dark halo of size ry radiative cools 
to form a disk of radial extent rj,. Before we consider the 
details of cooling in j |2.2| , we first need to estimate the 
collapse factor, r^/rd, which is based upon the conser- 
vation of angular momentum. The disk forms when the 
gas becomes rotationally supported. The protoquasar ac- 
quires its spin from tidal torques of its neighbors and N- 
body simulations (Barnes & Efsthaiou 1987) and analyt- 



(2) 



(4) 



ical studies (Heavens & Peacock 1988) indicate that the 
spin parameter of the virialized system, A^, to be in the 
range 0. 01-0.1. If the angular momentum is conserved 
and there is no exchange between the gas and dark com- 
ponents, the ratio of angular momentum in the gas to the 
halo would remain as Md/A4, so that 

Xv \EJ \MdJ 
where A = L\E\^/^ /GM^^^ is the spin parameter, L is 
the angular mometum, and E = —kiGM^/ry and Ed — 
-~k2GM^/rd are the binding energies in the halo and the 
disk respectively. If the halo has a constant density then 
ki = 0.3 and if it were a truncated isothermal sphere 
then ki — 0.5. Using the form for the circular velocity, 
Vc = y^GAi/ry, of a disk spinning in an isothermal halo 
and taking the angular momentum of the disk to be equal 
to sMdTdVc, where s is a geometrical factor (typically of 
order unity; s = 2 for an exponential disk) that depends 
on the mass distribution in the disk, the radial extent of 
the disk is given by 

Td = ry ( A_ I (5) 



Typically, the disk size is a tenth of the virial radius, and 
the collapse factor ranges from 10 to 20. For convenience 
we make the following definitions 

(6) 



Xd 



10 kpc 

rd 



1 kpc 

where rd is the radius of the disk. 



2.2. Cooling 



(7) 



In order that gaseous disks (with collapse factor estimated 
above) form in the halo where star formation and super- 
novae take place, it is important to examine whether the 
gas can be retained in the hosts in the first place. In this 
section, we examine the constraints on black hole hosts 
(those that can retain the gas) from star formation and gas 
loss due to supernovae. Consider a virialized halo which 
contracts due to cooling to a radius where it fragments 
due to self-gravity and star formation takes place. Below 
we list the conditions that specify when a halo condenses 
to form stars, whether the gas becomes unbound due to 
supernovae heating and finally whether the gas cools be- 
fore it can escape and hence remains trapped in the halo. 
The first two of the conditions were considered earlier and 
in this paper we introduce a third necessary condition pre- 
viously not considered. 

CO. Following earlier work (eg. White & Rees 1979, Rees 
& Ostriker 1977, Silk 1977), we state the condition that a 
luminous core can form in a halo 

tc < td (8) 
where tc is the cooling time and td is the dynamical time. 
The heating process by was examined in some detail by 
Dekel & Silk (1986, DS) in the context of dwarf galaxies; 
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they found that a condition of gas loss amounted to the 
virial velocity being below a certain critical velocity. At 
a time, i/, when the hot gas in supernovae shells signifi- 
cantly fills up the volume under consideration, the follow- 
ing constraints should be satisfied- 

Cl. As given by DS, the effective heat input by supernovae 
is 

E^tf) > \fgMvl (9) 

where fg is the gas fraction in the total mass and Vc is the 
circular (or virial) velocity in the halo. This implies that 
the supernovae heat input should be greater the binding 
energy for the gas to escape. 

C2. We find that an additional condition for gas loss is 
necessary, namely, that the time for the hot gas to cool 
should be longer than the escape time, tf. of the system 

> (10) 
where t,, is approximated by the time scale 

r" dr 

te^ I _ (11) 

■^0 ^J2{E - $) 

where E represents the mean kinetic energy per unit mass 
and $ is the mean gravitational potential per unit mass. 
In other words, even if supernovae heat input causes the 
gas to become unbound, it can still be trapped in the halo 
if it cools faster than the time required to escape, which 
is of the order of the dynamical time. 

In order to quantify these physical constraints, we 
have to calculate i?sn, the effective heat input by super- 
novae. The standard evolution of a supernova remnant 
goes through two phases- adiabatic and radiative. In the 
adiabatic phase (Sedov- Taylor), the radiation loss is negli- 
gible and the time at which the shock front radiates about 
three-quarters of its initial energy is given by DS as 
trad - 1.4 X 10^ yr, (12) 

where n is the ambient hydrogen number density in cm^"^. 
As a result, the input into the gas equals the initial en- 
ergy minus the radiative losses in the adiabatic phase and 
subsequently for later times, the gas is cooled by expan- 
sion. The cumulative energy input from the supernovae at 
a given time, will be dominated by those that have ex- 
ploded within a time trad before t. The star formation rate 
is taken to be a constant and is approximately, f*M/td, 
where /* is the baryonic mass fraction in stars and fur- 
ther, it is assumed that the star formation abruptly ends 
at t = td- The IMF of the solar neighborhood gives rise to 
one supernova per 200 Mq of newly formed stars. If ^200 is 
the number of supernova explosions per 200Mq of newly 
formed stars|^ with each explosion releasing eo = 10^^ ergs 

^ An analytic fit to the IMF (eqn(f .3.10), Shapiro & Teukol- 
sky, 1983) which is based on Bahcall & Soniera (1980), yields 
the ratio of the number of objects in the mass range 5 < 
M/Mq < 30 to the total mass in all the stars (with a cut- 
off at 0.1 Mq) to be 5 X 10"'^ per solar mass; this we have 
adopted as a fiducial value of the number of supernovae per 
solar mass of newly formed stars. 



of initial energy, then the energy input into the gas is then 
given by 

i?.„(0=eo/, 2o^-20o^/(t) (13) 

and the dynamical time for the system is a quarter of the 
oscillation period given by 

td^i v3 ^ ^-^ "^'^^ ^^^^ 

which is assumed to set the timescale for star formation. 
Also, f{t), is a function of time that is of order unity and 
its form is given by eqn (45) in DS 

{t/trad)[l - 0.14(t/t™d)17/5] t < trad 
0M + 0.58[it/tradf-^^ -I] t>trad 

The ratio trad/td 



fit) 



(15) 



0.01 y/h/OJ implies that /(<) is of 
order unity and < f{t) < 3.6 (for <t < td). The heat 
input from the supernovae that have exploded within a 
time, trad, of a given instant are the most effective in con- 
tributing to the heating. Now the shells of the supernovae 
will start filling up the volume. The time at which hot gas 
has a filling factor of order unity {t = tf), is estimated in 
the following manner. 

2.3. Filling factor of SNR 

We use the well-known simple expressions (eg. Spitzer 
1978) for the advance of the supernova through the Sedov- 
Taylor and snow-plow phases in a medium whose number 
density of hydrogen is typically n — (l/130)(A^io/a;^) = 
(1/130)(^)3, which is in the range of 0.01 to 1 H atom 
cm~'^. In the Sedov- Taylor phase, during which the to- 
tal energy in shock front is conserved, the advance of the 
shock front is given by 

r,t = 10^^ t^/^ n-^/^ cm, (16) 
where t is seconds. This phase ends when the temperature 
falls below lO^K and radiative losses are significant at a 
time trad, given by eqn ([l^). The radius at beginning of 
the radiative phase as given by this condition is then 

Trad = 40 pC (17) 

Next, the radiative (snow- plow) phase follows, in which 
the momentum is roughly conserved and the shock front 
advances according to rg oc i" '^^ (Chevalier 1974). The 
star formation is expected to occur when the gas in the 

halo shrinks to some size R { fbTy if it is isothermal or 

1/3 

^ fb ^^'^ cloud is uniform) and fragments due to 

self-gravity (Larson 1974). Taking the typical size of the 
remnant as given by eqn (|l7|), the total shell volume of 
the supernovae remnants in units of the volume occupied 
by the gas is 



Fit) 



t_ rl^^ ly^oo /.(Af/2OOM0) 
td 



10 



1 



i?3 

4/, -1/5 t 

- f* td 



(18) 



where eqn(|T^ and the density of the gas, Af(l — 
/*)/(47ri?'^/3) = 1.4 n nip was used. Since the star for- 
mation shuts off at t — td, the maximum value of the 
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ratio of the total shell volume to gas volume is F{t > 
to) = 10 (4/^/(1 - /*)) n-^l^. This implies that the 
filling factor of the supernovae shells in the galaxy will 
be weakly dependent on number density. One can take 
into account the porosity (which can be thought of as the 
complement of the probability that a given point in the 
volume is outside the N — 1^200 /* -^9 (5 x 10^) remnants 
of fractional volume, q = F(t) /iV, occupied by one shell) 
by(3 = l — (1 — (7)^, where Q is the filling factor, which is 
well approximated in the Poisson limit of large N by the 
formula, 

g(t) = l-exp(-F(i)). (19) 
It is clear that Q at large times is close to 1 (and nearly 
independent of density) and the hot gas fills the volume. In 
order to estimate the total energy input into the medium 
in eqn(p^), we calculate the time, t/, when Q ~ 1, -F(i/) = 
3 which leads to tj = F{tf)td/F{td) ~ 0.3^^ = 30 Uad- 

2.4- Gas loss criteria 

The solution leads to /(i/) ~ 3, from eqn([l5|). So the 
energy input into gas from supernovae works out to be 

Esnitf) = 3 X 10^5 ^ A ^200 Mio ergs. (20) 
The escape energy is given by 

1 „ . . 9 1 GfgM^ 

2 r,„ 



2 fa M vl = 



4 X 10 



55 (1 - /.) h Xfo 



0.8 0.1 



TTT ^^^^ 



(21) 



where 

Vc = 0.7 X 10^ ^/Mw/x^cui s^i (22) 
is the virial velocity. As a result the condition CI for gas 
removal can be written as 



sn 
.,2 



3/, 



1^200 > 1 



(23) 



1 - /* Mw 

where 2E,n/{fgM) where fg = 1 - /6(1 - /*) is 

the gas fraction in the halo mass not converted to stars. 
It is clear that enhancing the star fraction increases the 
energy input into a smaller gas fraction and hence vl^^ is a 
monotonically increasing function of /*. Also, halos of the 
same mass had deeper potential wells in the past which 
trap the gas better. Similarly, more massive halos clearly 
have deeper wells and gas loss is less likely as seen from 
eqn(H). 

In appendix ^ we calculate the difference in gravita- 
tional binding energies in the initial configuration of an 
isothermal halo and a final one consisting of an exponen- 
tial gaseous disk in a halo consisting of stars and dark 
matter. The gas taken to be roughly near the virial tem- 
perature, 

T. = ^^lO^K (24) 

3 kit Ty Xy 

where the mean molecular weight, fi — 0.6. The cooling 
proceeds through line and free-free emission and we take 
a cooling function, A(T) — 10"^^ ergs s~^ cm~^A23(r), 
provided by Sutherland & Dopita (1993) for a metallicity 



of [Fc/H] = —4. The cooling rate in the virialized halo 
during the contraction is roughly given by 

Eh^An^Ne^ 4.5 x lO-^^ A23(r„) M| x^^ ergs s-\{25) 

where is the number of electrons. Now, we can estimate 
the cooling time taken for the system to cool and fragment, 
given that the source of thermal energy in the gas is one- 
half the change in gravitational potential energy from the 
virial theorem. The cooling time, tc, is given by 
_ AI^/2 



(26) 

Eh{Ty) 

where Ty — lO^Mio/xy and plugging in the expressions 
for AW from appendix |^ the condition of core conden- 
sation CO, thus, can be expressed as 



/0.08A 


y/MioXy 


'0.3 1" 


V /. ) 


94 


Xy 2 



A23(T„)"1 < 1. 



(27) 



After fragmentation and star formation, the gas is 
heated up and the hot gas cools in a time, tg = 
Esn/E{Tsn) where Ts„ = v'^„iJ,mp/{Skb) = STy = 
10^ S Mio/xy is the temperature of the heated gas and 
^ = '"'in/'^c = {xyiy2oo/Mio) 3/,/(l - /,). Now, we need 
to compare this to the escape time, as defined in eqn 
(pi]), when the condition CI is satisfied or a quarter of 
the oscillation period when the system is bound. Now, 
E = Esn/ifgM) - vll2 = vl{b - l)/2, is the mean en- 
ergy of a gas particle in the system and is zero when the 
escape velocity equals the critical velocity. Using the po- 
tential for an isothermal sphere that is truncated at r^,, 
$(r < r„) = ■i;2(ln(r/r„) - 1), and Vyjvc = {2/TT)td, we 
obtain 

~ 1 + .5, 



te{5) = td\l -CXp{ 



1-Erf 




5 > 1(28) 



which is valid for 6 > 1, for which the supernova heat 
input, Esn, is such that the resulting mean energy of the 
particles is positive. These particles can escape only if the 
cooling time is longer than the escape time. The escape 
time, te{S), is a slowly decreasing function and is of order 
td (See Fig I). For (5 < 1, the gas remains bound and only 
the condition of core condensation, CO, applies. Hence, the 
condition C2, tg/te > 1, can be expressed as 



where the effective dynamical timescale at (5 = 1 is taken 
to be td and te(l) — 0.418 td, so that the curve C2 is 
continuous with CO. This is seen in the cooling diagram, 
presented in Fig || where the number density of hydrogen 
in the halo, n(cm^'^) = Mio/ {150x1), is plotted against 
the virial temperature, Ty and shows the curves given by 
CO, the lower curve, CI, the vertical line, and C2 which 
is the upper curve in the region to the left of CI. The 
cosmological parameters were set to /170 = 1,^ = 1, 
and fb = 0.1 and the halo parameters were chosen to 
be Xy = 0.075, /» = 0.2, and 1^200 — 4. The shaded region 
indicates the halos that have collapsed but star formation 
has induced gas loss and they ultimately resulted in gas 



A. Mangalam: Formation of Proto-quasar 



7 



12 ■ 



2 4 6 8 

SN ENERGY INPUT/HALO BINDING ENERGY (6) 



Fig. 2. The dependence of the escape time of the hot gas, 
in units of the dynamical time in the halo which is defined 
as one-quarter oscillation period, on S — vln/v'^, the su- 
pernovae energy input in units of the halo binding energy. 
A factor of 10 in the heat input reduces the escape time, 
tf., by factor of 2. 




4.25 



4.75 5 5.25 
Log Tv (K) 



Fig. 3. The shaded region bounded by the constraints, CO, 
CI, and C2 in the cooling diagram contains the halos that 
have cooled and contracted to a radius of fragmentation; 
however, the supernovae heated gas in these systems have 
escape times shorter than the cooling time, resulting in 
gas-poor galaxies. The parameters chosen here are /, = 
0.2, 1^200 = 4, and A^, = 0.075. The halos above C2 for 
5 > 1 (on the left of CI) and above CO for 5 < 1 can 
trap the gas and hence are the candidate hosts of central 
massive objects. 



poor systems. The halos above this region retain the gas 
that could form the massive central black holes. 

The corresponding range in mass for a given redshift is 
shown in Fig ^, for the choices of the parameters, z/200 — 
4, the supernova efficiency, and /» = 0.2, the fraction of 
baryons in stars. Clearly for large redshift, the mass range 
increases and this is due to the deeper potential well which 
can retain the gas better. There is a crucial question of 
whether all the gas is retained or if some fraction is lost 
in a wind. A supernova efficiency in the halo of 
4/* 



^■200 



■> m^Mw/xl = 4.3(1 + zcf/^, 



(30) 







11 - 



cn 
cn 

o 

< 9 - 



O 




1 2 
COLLAPSE 



REDSHIFT Zp 



Fig. 4. The halos in the shaded region in the — Zc 
space can trap the gas; the parameters chosen are /* = 
0.3, J^200 = 4, and A„ — 0.07. 



from eqn(^9|), would be required to drive winds that would 
result in mass loss and this is considerably more than the 
corresponding estimates in our galaxy, for 2. 

If the stars form after the disk forms, the cooling rate 
in the disk would be given by 



i?. = 3x10^4 A,, f 



M 



...?o^."'(W*)'gergss-i(31) 

and the cooling time would reduce by a factor of a thou- 
sand corresponding to enhanced density. Although only 
the cosmological abundance of [Fe/H] = —4 was used, 
the supernovae explosions will also enhance line cooling 
from metals injected into the medium; this would reduce 
the cooling time estimated here. 

3. Magnetization of the disk 

In the large number of the halos which trap the gas as 
given by the conditions C0-C2, a gaseous disk forms with 
a radial extent given in 



2.1. Further, it was seen in S2.E 



that the supernovae shells fill the core volume at the time 
of star formation, and the small scale magnetic fields are 
dragged with the gas as it settles into a the disk. Now, we 
consider the question of whether the field strength is large 
enough to provide a significant viscous stress. We take 
the gas to be initially dominated by the gravity of the 
dark halo and assume the following logarithmic potential 
(Binney & Tremaine 1987, eqn (2-54)) that obeys the flat 
rotation curve 



(i>(r, z) 



(32) 

where tq is the radius of a compact region, < 1 and 
Vc is defined in eqn (|2^). If the vertical equilibrium was 
a result of balance between vertical gradient in the total 
pressure, P = C^^P (which represents a sum of magnetic 



8 



A. Mangalam: Formation of Proto-quasar 



and gas pressures and the density scale height is H), and 
—dz^, we obtain 



H 



Assuming a dominant isothermal halo, we obtain the half- 
flaring angle of the disk, 

(34) 



' d 



where Q — (1/8) (ql/O.S), and the speed of sound 



was taken to be Cc 



firrip 



1.6 X 10^ 



cm s 



for 



a temperature of lO^K, where the cooling curve drops 
significantly. 

The supernovae inject the medium with magnetic flux 
and next we estimate the typical strength of the small 
scale field. The calculation in §2.3 shows that the filling 



factor of the shells in the total gas volume is nearly unity 
for a wide range in density and therefore the supernovae 
hot gas fills the medium and hence quasi-uniformly mag- 
netizes it. The volume of the flared disk works out to be 



Att 



r3 H 



= 1.67 X 10® pc^ 



(8 9) 



4.6 



0.075 J V1 + --, 
The Crab Nebula of size 0.8 pc (or volume 2 pc'^) has 
3 X 10~^G fields (see §0 for a discussion of supernova field 
strengths). By freezing the flux in a Crab volume of 2 
pc'^ to a volume occupied by one shell, the magnetic field 
strength in the disk for N = 10SA^ioi^2oo(/*/0.2)(/b/0.1) 
remnants is 

fVd/N\-'/' 



B,„ = 3 X 10"* G 



1.7 X lO^^G 



1 + Zc 



4.6 



2/3 



A 

0.2 



2/3 



A 

0.1 



2/3 



0.075 



(36) 



By choosing a typical set of values (A„ — 0.05, 1^200 ~ 
4, = 0.5, 9 = 1/8) the field strength turns out to be as 
high as Bsn = 170/iG. 

The mean number density of hydrogen atoms in the 

flared gaseous disk, for /i = 0.6, is given by 

rd\ -3 

h) 

7/2 



rid 



7.5 (1 - U) Ah -TT 



fl + z, 



M 



2/3 
10 



A 
0.1 



(37) 



= 150 cm 

0.8 V 4.6 

where (||) and (|) were used and the values of A^, — 0.075 
and 9 = 1/8 were taken. Outfiows from O and B stars 
could also magnetize the gas, and by using flux freezing we 
estimate the field strength due to winds to be (Bisnovatyi- 
Kogan et al 1973; Ruzmaikin, Shukurov & Sokoloff 1988) 

= {p/p^f'^Bs = 5.2^G (^)'^' • (38) 

where the estimates of the density at the base of the wind, 
~ 2 X 10"^'^ g cm~^ and the field at the surface of the 
star, Bg ~ 4G are used. The field expelled by massive stars 



are unlikely to pervade the volume, and the major contri- 
bution would be from supernovae. Although the number 
of massive stars are of the same order as the number of 
supernovae using the Salpeter IMF (estimated by inte- 
grating the mass range above IOM0; see the footnote in 
§2.2), the smaller fluxes of the wind (^ BsR^ where i?* is 



the radius of the star) render the effective field strength 
to be weak. The key point here, is that for typical values 
of the halo parameters, Bsn is about 10"** — 10~^ G which 
is a factor of 10-100 higher than the value in our galaxy 
(Vd = 6 X 10"pc3,iV - 3 X 10^, Bsn = 3AiG; due to a 
smaller value ofV^/N of the supernovae shells. 



4. Accretion to a compact region 

Here we consider different viscosity prescriptions, namely 
direct magnetic stress, the phenomenological a viscosity 
(Shakura & Sunyaev 1973) derived from magnetic fields, 
and angular momentum transfer mediated by self-gravity 
induced instabilities, and calculate the corresponding ac- 
cretion timescales. The time dependent disk accretion is 
described by the conservation of mass, radial momentum 
(where all other forces except gravity are neglected, and 
Vr/vrj, <C 1) and angular momentum 

rdtY. = ~dr{rY. Vr) (39) 
LJ^ = (l/r)9^$ (40) 

Y.VrdAr'^oj) = l/r drir'^Wr^,), (41) 
where S(r, t) is the surface density, $ is the gravitational 
potential and 11^0 is the vertically integrated stress. We 
take ro as the inner radius of the disk fiow and the outer 
edge of the compact region while 



uj(r) = - , 

By integrating eqn (^ij), we obtain for steady flow 



Muj 1 



2r^uj 



(42) 



(43) 



The viability of the various alternatives for the viscous 
stress can be assessed by estimating the accretion rates. 
In the initial phase of accretion the potential is dominated 
by the halo as calculated in §|[ while in the later phase, 
self-gravity dominates and the resulting flow is calculated 
in §^. 

4.1. Magnetic Stress 

We consider now the form of the stress tensor that is en- 
tirely due to magnetic fields. We assume that the magnetic 
field injected is largely in the form of small scale loops 
that are frozen into the plasma. Further, we expect that 
the processes of compression and advection preserves the 
form given by 

B = Bs{p/psf'' (44) 
that is likely to be valid for small scale fields in the non- 
dissipative limit. The Lorentz force on the plasma is 



(J X B)^ = 



1 

47r 



1 



■BrdrirB^) + B,d,B^ 



(45) 
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where the first term on the right hand side is the lo- 
cal shear stress and the second term is negligible in the 
initial phase (the vertical average of small scale fields, 
Bz{dzB^) w 0), but it could be important during a later 
phase of accretion where magnetic braking can operate via 
built up, large scale Bz, by angular momentum transfer 
to the parts external to the compact region. Here, we take 
only the local shear stress « B^H/{2-k) and this was as- 
sumed to be initially at sub-equipartition levels as a ther- 
mal pressure (C = 1) was used for calculating the half- 
thickness of the initial disk in eqn (^4|). One can verify 
this through a consistency check by using the definition, 
C = 1 + (STTpCg), which determines the half-thickness, 
H oc from the condition of vertical equilibrium, (34), 
Bsn C^^^'^j from eqn (^6|), and p (x ^^^/^ from eqn (37). 
It follows that 

-1/6 



0.03 C 



i = C 



(46) 



and for a reasonable choice of parameters (1^200 = Ij /* = 
0.2, Zc = 3.6, A„ = 0.075) this results in C = 1.03. Low 
values of 1//3 = Bs„/(87rpc^) ~ 0.03, are typical for the 
range of interest in the parameter space. 

Next, we estimate the accretion time scale and hence 
the viability of magnetic accretion in the steady limit. 
For typical values of the parameters assumed here (A^, = 
0.075, 9 — 1/8) the accretion rate turns out to be 



M 



27rH 



r0 



0.2 Bi 



B^H 



4.6 

1 + z, 



-) 



]B^ 



M, 



yr 



(47) 



using eqns (^, |2^, ^) . Taking the initial field to be from 
supernovae expulsions (B = Bsn) from eqn(p6|), this im- 
plies a time scale of magnetic accretion. 



M 

2 X 10^ yrs 



4.6 



M 



-2/3 
10 



-4/3 
^200 



ni;, 



o.i 



4/3 



(48) 



1 + Zc 

O^x 4/3 
/* / 

Having obtained this fast timescale to accrete lO^Af© of 
gas, we proceed to calculate the detailed form of HJ.*^. The 
thickness of the disk is given by eqn (^^, the balance 
of the total pressure gradient and vertical gradient of the 
background potential. Combining eqns ( ^ , (^) , and p ~ 
E/(2iJ), we obtain 

H^, = -(S?/4.)p;4/3e4/V-V30-i/3 (49) 

where Oq = 29{Aiio/xvy^^ is the flaring angle of the full 
thickness of the disk. Initially, the magnetic pressure is 
lower than the thermal pressure but as the matter sinks 
into a compact region the magnetic pressure is expected 
to dominate the vertical pressure gradient. 



4-2. a viscosity 

If a small scale dynamo operates quickly (Kasantsev 1967, 
Kulsrud & Anderson 1992) then it will saturate near 
equipartition values (as is well known from simulations 
- Hawley, Gammie &; Balbus 1995 and references therein) 
and an appropriate form of the stress can be described 
in terms of a prescription of the form H^^ — —2HassP, 
where P = Cpc^ is the total pressure, which is propor- 
tional to the gas pressure at equipartition. The accretion 
time scale is expected to be similar to the one obtained 
earlier. The dependence on the halo parameters can be 
expressed in the isothermal limit as 



„ ^ o 0.1 / 4.6 

- 2 X 10« yrs — — 

Uss \1 + Zr 



1/2 



0.1 



X M 



-2/3 
10 



0.^ 



h8- 



(50) 



,1-/*, 

Note that, although we use a direct magnetic stress in our 
calculations, we record for comparison, the detailed form 
of the a prescription in Appendix 

4-. 3. Gravitational instabilities 

Cold, thin rotating discs are known to be unstable and the 
basic stability criteria was provided by Toomre (1964). 
For a uniformly rotating isothermal disk (Goldreich & 
Lynden-Bell 1965) the criteria for local stability is given 

by 

Qt - g| > 1.06. (51) 
We find that 



Q 



TTvlH 



M 



-1/3 
10 



4.6 



1/2 



0.8 



0.1 



(52) 



GMd 8 \l + zcj 1-/, 

where the values {0 — 1/8, A« = 0.075) were assumed. 
So the disk is unstable to gravitational instabilities and 
it is possible for angular momentum transport to occur 
through this process. Lin & Pringle (1987) estimate an 
effective kinematic viscosity for gravitational instability 
from — LuP, where the critical shearing length, i = 
GEcj^^ was taken to be the maximum possible size for the 
instability. Here we make a more conservative estimate by 
introducing the parameter, 
given by 

H^0 — ag ^ rdrUJ = —ag T,^ lo^"^ if Qt 
The corresponding timescale of accretion is 



J, < Ug < 1, into the stress 



< 1. (53) 



lO'* M, 



© 



2tt H^, 

r<p 



(54) 



lO^yrsTW^/^ 



'10 



4.6 

l + Zc 



1-/* 



0.1 

X 



0.1 



■ m 



5. Self-similar evolution of the disk in a 
background potential 

Having demonstrated that the time scales of accretion are 
quite fast, we now proceed to calculate a detailed model of 
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self-similar evolution of a disk from the diffusion equation 
obtained from the conservation laws ( p9[ - |4l| ) 

9tE = -dr ( . dr[r^Ur4,]) , (55) 

r \dr{r''u;) ) 
where the viscous stress can be parameterized as 11^0 — 
—K-iI^t'^ and in addition the rotation law is assumed to 
be of the form lo = Kxr"^ . This very useful formulation of a 
self-similar form is due to Pringle (1981) but only partic- 
ular analytic solutions to the diffusion equation has been 
reported for the specific cases of (a = —3/2,6 = c = 3; 
Lin & Pringle 1987) in the context of accretion of a pro- 
tostellar disk onto a point mass via gravitational instabil- 
ities and (a = —3/2,6 = 5/3, c = —1/2; Cannizzo, Lee 
& Goodman 1990 (CLG), see Appendix ^ of this paper) 
in the context of disk accretion of a tidally disrupted star 
onto a massive black hole. Note that in CLG, the scaling 
law for the viscous stress by the closure of the conditions 
of local dissipation in an alpha disk in a Kepler poten- 
tial and vertical equilibrium. Here, an analytic solution to 
the general problem of the type (n^^ cx S'' r'^, a; oc r°) is 
presented so that possible viscosity mechanisms discussed 
earlier and expressible in this way, can be explored within 
the same formulation. In the magnetic case given below, 
the viscosity scaling is due a magnetic stress, the flux- 
freezing condition and vertical equilibrium in a cold disk 
in the background halo potential. A solution for an alpha 
disk with local dissipation with a general rotation law is 
provided in Appendix The general solution presented 
below has a larger utility in contexts other than one con- 
sidered here. 

If 6 = 1, the equation is linear and the general solution 
is easily found. Proceeding generally, under the assump- 
tions of self-similarity for (6 ^ 1), one may write the the 
surface density in the following form 
S = Eo(Vio)''ff(0 rj=r,{tlt^r, (56) 
where ^ = r/r/, and is the associated radius scale. We 
set the constants ^ 

to = ^rl+^-^Y.^, Ma = 2TTlEorl I = [ ^gm^ (57) 

J^2 Jo 



where Md is the initial disk mass. Here we seek a particular 
solution when there is no external torque, which implies 
the total angular momentum, J, of the disk is a constant. 
Using the scaling relations above that are implicit in eqn 
( |55| ) and 

J = 2^ifirf+"S]oj, J - / egm^ (58) 



it follows that a = {Ab + ab - 2 - c)~^ and /3 = -{a + A)a. 
At this point we note that the disk edge travels outward if 
2 + c < 6(4 -I- a). Substituting into the form for the surface 
density, as given in (56), and simplifying (p5[), we obtain 
the following ordinary differential equation, 
1 



a^g' + Pg 



rd, 



(59) 



(a + 2)e ^ 
After some algebraic transformations, one can integrate it 
once to obtain 

a{a + 2)gr+' + (a - c)e+'g'' = be+'g'-^g' + c,. (60) 



Now, we apply the boundary condition that the density 
vanishes at the disk edge, ie., 

^{r = rf)=Q; .9(1) = 0. (61) 
Moreover, if 6 > 1, which is the case for the examples 
considered here, then c\ — 0. By rearranging terms and 
integrating, we obtain the following solution 

g{i)-&-'^^"{i-e^'^y~'^ ■ (62) 

The time constant can be evaluated with ~ rd to be 



Md 



6-1 



Md 



m 



hS 



(63) 



|Hr0(So,'"d)| 27r/ 
where ta is the accretion time scale calculated in §4.1 and 



§4.3 in the steady case using Hr0(Md/(27rr^), r^). The rate 
at which mass sinks into the center is given by 

{a+2)a' 

1 



M,it) = Md 



(64) 



The accretion time scale which is the time that transpires 
when a fraction e = THhs/Mds of the disk mass falls in is 
given by 

Ta = to ((l-e)(^ - I 



tn 



' (< 



1 



(65) 

where ta represents the timescales tm , tg or ta , and indi- 
cates that the estimates derived earlier are modified by 
geometric factors. Now we consider the particular case of 
magnetic accretion (a — —1,6 = 4/3, c — — 1/3, e = 1/8, 
see §4.1 where the value of e is chosen for a typical case 
where the disk mass, Mds — 8, and THhs = 1) which has 
the solution 

a = 3/7, = -9/7, T,„ - 1.01 

3 



gmiO = f 



-1/2 



3/2 



(66) 



Similarly, the accretion due to gravitational instabilities 
(a = —1, 6 = 3, c = 2, e = 1/8, Mds = 8, and Hl/ig = 1, see 



§4.3) has the solution 

a = 1/5, /3 = -3/5, Tg = 0.54 ig, 

ff.(C)=rMl-?]'^'- (67) 
The disk structure of these solutions are shown in Fig. |^. 

The self-similar solutions of this kind to (|5^) are known 
to develop at large times in numerical simulations with a 
variety of initial conditions (Lin & Pringle 1981). 

Now we determine the regime in parameter space 
where the halo dominated flow can occur. This is given 
by the condition that a,t t = to and r = rd, the halo dom- 
inates self-gravity or 







d£d 




> 




rd 


dr 


Mh 


> 


Md 


Tdrv 


Ir'd 



Va^Ji(an) / Jo(^an)ff(OfdC, 
"'0 



(68) 



where is the circular velocity due to the halo alone, 
the disk potential was expressed as a Bessel transform of 
the surface density and a„ are the zeros of Jq. The col- 
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Fig. 5. The structure g{^) of the self-similar disk, S — 
^0 (^)^ 5(0 where ^ = r-/r/,r/ = r, (^) . The mag- 
netic solution, (?m(0 ~ ^^^/^ [l^C''^^]'^ is shown by a 
solidline, and the solution of the disk with gravitational 
viscosity, gg{S,) = S,^^ [1 — C]^^^? is shown by a dashed line. 



lapse factor r„/rd = sVfci/Ai, = {j/l)vki/Xv from § |2.1 
Mfi/Md = (1 — /g)//^ by definition, leading to 
A.(l-/<,) 



> 



(ean)5(0^d^ (69) 



For the solution (|66|), the magnetic case, the RHS of the 
above equation works out to be 0.44 and for the gravita- 
tional instability case, the solution (67), the corresponding 
value is 0.46. Hence this condition can be written as 

A.(l-/9) 



> 0.32 



(70) 



where ki = 0.5 for a truncated isothermal sphere was 
used. For a reasonable range in the parameters (A- 



0.05 - 0.1, /g = 0.05 - 0.08) the above condition holds 
good (Au(l — fg)/fg is in the range 1-4); as a result, the 
initial accretion flow is expected to be halo dominated. As 
the mass accretes into the center, the spin deviates from 
w oc 1/r and gradually increases. The self-similar solutions 
are valid only up to a point beyond which the self-gravity 
due the disk and the central mass dominates the poten- 
tial. The time of transition to a self-gravitating flow can 
be estimated by seeking that 

dr r 

Now expressing time in terms of the e = Mc/M^ in eqn 
(p^), we obtain after some algebra 

'a^Ji(a„(l - e)y) / Jo(Can)5(OCdC 



> 



A.(l-/s) 



(72) 



W - e) " faVk~i ' 
where y — r/r^. Clearly, this condition docs not hold close 
to a compact central region (defined by r < tq) where it 
is dominated by the central mass. The solutions consid- 
ered here are a good approximation for the region be- 
yond y = 0.1 for e < 0.3 and provide a sufficiently accu- 
rate description (See Fig. At large times, a Keplerian 
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Fig. 6. A density plot of the ratio of self-gravity 
(disk and the central mass) to halo gravity defined as 
j_ _|_ GM^{t) ^ ^ giiowing the evolution of a mag- 

netized disk in the halo potential as the central mass 
increases to e = M^/M^ = 0.3. The horizontal axis is 
the in units of the initial disk radius, r^. The value of 
A^(l — fg)/{fg\/ki) was chosen to be 3; the halo gravity 
(in units of v^) dominates at radii outside a central region 
of 0.1. 

flow into the compact region can be assumed to occur 
and one can use a self-similar flow again with a — —3/2 
and the corresponding magnetic stress taking into account 



Vc oc r 



P 



a+l 



in eqn (p4) (in combination with (M 



= S/(2iJ)) leads to H™, cx S^/a^a/s ^ ^^/^r 



, and 
Sim- 

ilarly H^^ (x S^r"^" = S^r^ from eqn (|5^). However, the 
estimate of the accretion timescale is not expected to be 
very different from the one derived earlier. The key result 
is that about a fraction, 0.3, of the disk mass can be trans- 
ported into a central region, which is a fraction, y = 0.1, 
of the initial disk radius within the time given by the so- 
lutions (|^. For typical values, this implies that about 
2 X lO^M© of gas sinks into a t-q = O.lrd ~ 100 pc region 
in a halo (with Alio = 1,A^, = 0.075) in a timescale of 
2 X 10^ yrs. 



6. Disk evolution in a self-gravitating regime 

We now examine the evolution of the disk in a gravity 
field that is entirely due to itself. This flow can occur after 
sufficient accretion of mass into a compact region of radius 
ro; or, alternatively at the time of the formation of the disk 
if the eqn ^ is not satisfied and \v{l-fg)/fg< 0.32. The 
problem of self-gravitating accretion flow is complicated 
by the coupling of Poisson's equation to the momentum 
and continuity equations. Clearly, its evolution has to be 
treated differently from the preceding case of a prescribed 
background potential. A detailed model is deferred to a 
work in preparation (Mangalam 2001); below, we consider 
a useful simplified version of the problem by assuming 
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a particular form of the density distribution (see Field 
1994). 

We assume a Mestel (1963) disk where the self- 
consistent density distribution with potential is entirely 
due to self-gravity, is of the form 

E(M) = 5*ii, (73) 

where the time dependence appears only in the rota- 
tional velocity. This is a an interesting disk which is self- 
consistent, taking into account the most relevant physics. 
Taking 

V4, = vo x{t) and r = tq Xi {t) x (74) 
where Vq — GMc/tq, and Mc is the mass out to tq. We 
see that by assuming a self-similar evolution of the disk, 
the mass out to a given x should be independent of t and 
hence it follows that 



Xi = X ^ and E = S, 



X 



(75) 



where Ec = Mc/(27rrQ). From the continuity equation. 







(76) 



rdtTi + dr{rT,Vr 
we find 

Vr = -Ax'^xrox. (77) 
Substituting this and the self-similar forms given above 
into the angular momentum equation, 
mv^ + ^{vr/r)dr{rv^) = {l/r^)dr{r^nr^) 

-Svo^cX^X = --^{x^Ur^), (78) 
X ax 

we obtain 
3 

= ~-^voT,crox^X, (79) 
which is independent of x. So far no specific viscosity mech- 
anism has been invoked- the form of II^^ above is neces- 
sitated by the prescription of a Mestel disk. If a magnetic 
stress is assumed and H^^ = — kB'^ H / {2Tr) , where /t is a 
factor of order unity and B is given by the vertical bal- 
ance of magnetic pressure and gravity, B — 2'Ky/G'E = 
2-K\fGT,cX^ I X. It follows that the half thickness, H (x x^, 
and can be expressed a,s H — Hqx^x^- Furthermore, from 
the flux-freezing condition, B oc [r^iJ]^^'^, , it is seen that 
p = —2. Equating the form of n^.^ from eqn(|79|) to a 
magnetic stress, kB'^H/{2-k), and writing (after taking 
e = 1/8, A^, = 0.075) 

1 Vf)rl ro 1 ro 



2GMaH Vo 2k Ho 
6 / ^o/rd 



10 



4.8 



0.1 J ' 1 + z. 

we see that 
3t" X X = X 

which leads to the solution 
r 



X(i) = 



t 



1/3 



(80) 



(81) 



(82) 



where x(0) = 1 was taken as the initial condition. So 
the disk spins up rapidly and shrinks to a smaller radius. 
Clearly, the solution is no longer valid when it is rela- 
tivistic. The self-similar collapse of the compact region is 
only a sketch but nevertheless the collapse timescale, r. 



suggests that formation of a black hole is extremely rapid 
(- 10<^ yrs). 

7. Summary of the results and discussion 

We summarize our results as follows: 



1, 



We considered star formation and supernovae feedback 
on the remaining gas in the halo based on the frame- 
work of DS. We include a new and necessary condition 
C2, namely that escape time for the hot gas be shorter 
than the cooling time and find that the condition for 
gas loss can be expressed by eqns ( p7| , p3| p9| ) which is 
depicted in Fig ^. For a typical choice of and ly, Fig 
^ shows the allowed range for Halo mass in terms of 
redshift. There is a sharp decline in the allowed range 
beyond collapse redshifts of = 2, as the potential 
wells formed at earlier epochs are deeper and trap the 
gas better. 



2. It was seen in §2.3 that the hot gas from supernovae 



has a filling factor of nearly unity and hence the small- 
scale magnetic fields occupy the disk. It was shown in 
§H the field strength is significant (10 — lOO^G), based 
on typical values of supernova efficiency, v, and star 
fraction, 

In we examined magnetic stress and viscosity due 
to self-gravity induced instabilities that would operate 
in a collapsed disk. The accretion timescales were esti- 
mated for a direct magnetic stress, eqn( [l8|) , an a pre- 
scription, eqn (]50|), and self-gravity, eqn(|55|), in terms 
of halo parameters {M, Zc) and star formation param- 
eters (/*, ly). The timescales are all short compared to 
the cosmic time. 

A general solution for a self-similar and time depen- 
dent accretion flow for a viscous stress of the form 
11^0 oc Ti^r'^ in a prescribed potential, ui (x r"', was ob- 
tained in §||. This was applied to the specific cases of 
magnetic accretion, eqn(|6^) and and gravitational in- 
stabilities, eqn( |67|) . The structure of the resulting disk 
is shown in Fig. |5| and the timescales are within geo- 
metric factors of their estimates. The disk eventually 
becomes fully Keplerian. The condition of halo domi- 
nated flow is given by eqn(fO|) and the transition to a 



self-gravitating flow is given by eqn(|72[). This solution 
is valid for dominant halos in the initial stages for the 
outer parts, y > 0.1, oustide a Keplerian compact re- 
gion for up to time when a fraction of the disk mass, 
e- 0.3, falls in (see Fig|). 
5. A self-gravitating Mestel disk that is spinning up as it 
is collapsing self-similarly, is described in We apply 
this to a compact region where the pressure is assumed 
to be due to magnetic fields. The time scale of collapse 
turns out to be the rotation time of the outer radius, 
which is a few million years. 

We now discuss some of the issues involved: 
Supernovae feedback 
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In a pioneering work, DS, showed that with a model 
of protogalactic gas in a halo reproduces the observed re- 
lations very well with an initial CDM model. Further, the 
condition for gas loss was given by the condition that the 
energy input from supernovae is more than the binding 
energy, CI and the cooling condition for gas contraction 
and star formation. In this work, we specify an additional 
necessary condition, C2, that the escape time be shorter 
than the cooling time. As indicated in Fig. ^, this restricts 
the gas loss to the halos in the shaded region. The exact 
shape and location of the region is subject to the choice of 
parameters (i^, /*), but the key point made here is that the 
entire region to the left of CI cannot be considered as a 
gas loss zone, even if heating by background UV photons 
were included. It is worth investigating this question in 
more detail by simulations that include the hydrodynam- 
ics and cooling of the hot gas. Based on this study, one can 
conclude that most of the dwarf galaxies involved efficient 
conversion of stars or high efficiency of supernovae, since 
that would place curve C2, higher on the cooling diagram 
accounting for the prediction of the location of the dwarf 
galaxies made by DS. The constraints imposed here, how- 
ever, account for the presence of black hole systems on the 
left of CI (and above C2). It also explains a sharp decline 
in quasars at epochs later than z = 2. 
Magnetic field strength in supernovae shells 

In estimating the field strength B^n in all the shells 
were assumed to have 10"'' G when they are young (size 
of 0.8 pc) and flux freezing was used to estimate the values 
at a given size. This field has been measured in the Crab 
nebula and in other young supernova remnants like Tycho, 
Kepler and Cas A, the field strength in the shell is inferred 
to be 10-3 - 10-4 G (Strom & Duin 1973, Henbest 1980, 
Anderson et al. 1991). Simulations by Jun &: Norman 
(1996), show local turbulent amplification by Rayleigh- 
Taylor and Helmholtz instabilities but claim overall sub- 
equipartition strengths. An analysis of Cas A X-ray and 
radio surface brightness at high resolution, show that there 
is a strong correlation, implying a possible equipartition 
between the field (^ mG) and the hot gas (Keohane, Got- 
thelf, & Petre 1998). The field amplification can be from 
other sources as well. For example, the magnetic field in 
the Crab has been wound up in the body of the nebula by 
the pulsar and the magnetic energy is a few a percent of 
the spin energy (Rees & Gunn 1974). High resolution sim- 
ulations of a supernova including rotation and expanding 
in an unmagnetized medium, which to our knowledge are 
unavailable, are best suited to the answer the question of 
the initial field strength in the disk. The small scale fields 
that are injected in the medium could be amplified to sat- 
uration levels by a dynamo operating in the disk aided 
by differential rotation and turbulent motion. Since this 
amplification was not taken into account (no assumption 
other than flux- freezing was made for the field evolution), 
the field strengths derived from supernovae which is in the 
range 10-100 fiG, could be an underestimate. 



Large scale fields 

Strong fields of a few are estimated from Faraday 
rotation of background QSOs by damped Lya absorptions 
systems, which are thought to be proto-galactic disks, as 
early as z 2, imposing constraints on the kinematic 
dynamo (Kronberg, Perry, & Zukowski 1992). A primor- 
dial origin for the large scale field in galaxies has been 
speculated upon (Rees 1994, Ratra 1992 and references 
therein), but the field generated by these mechanisms is 
much smaller than the magnitude required for a prime- 
veal hypothesis. Large-scale magnetic fields in the galaxy 
are thought to be generated by a turbulent dynamo from 
a weak seed field 10-^^G). But, it has been argued 
that magnetic energy at small (eddy) scales builds up 
much faster (at eddy turnover rates) than the mean field, 
so that the kinematic dynamo shuts off before the large 
scale fields amplify to observed levels (Kulsrud & Ander- 
son 1992, Kasantsev 1967) if the initial seed is weaker 
than 10-^G. Subramanian (1998) argues, however, that 
the fields are intermittent and do not fill the volume. The 
small scale flux ropes are at saturated values but the mean 
energy density is at sub-equipartition levels allowing the 
large scale field to grow to equipartition strength, aided 
by ambipolar drift. This question is still unresolved, but if 
the supernovae field can provide a fairly strong seed field 
(> 10-^G), then a kinematic dynamo can operate to am- 
plify the fields to fiG strengths in a galactic disk. A sphere 
of radius L, encloses (L/rgn)^ remnants where r^n is the 
size of a remnant. The rms value of the field at a scale 
L > rsn would be roughly ^ i?s„(rs„/L)'^/^. At the scale 
of the compact region, L = rQ — lOOpc, the large scale 
seed field is of order 10-^ G. 
Angular momentum transport 

There is an angular momentum transport mechanism 
required to form a seed mass in a central region and 
one for the compact central mass to further collapse; the 
two schemes, in principle, can be different. Loeb & Ra- 
sio (1994) concluded from hydrodynamic simulations that 
fragmentation occurs halting a direct collapse to relativis- 
tic scales and suggested that low spin systems could form 
supermassive seeds of disk or star geometry which would 
contract under radiative viscosity. The picture proposed 
here is that direct magnetic stress or gravitational insta- 
bilities in the disk can transport matter into a compact 
region on a dynamical timescale. The scheme of gravita- 
tional instabilities hinges on Qt < 1 in the disk whereas 
the direct magnetic stress depends on the supernova effi- 
ciency and field expulsion. The presence of high magnetic 
energy density from MHD turbulence has been invoked to 
explain high velocity dispersions of H I clouds in galactic 
disks (Sellwood & Balbus 1999). 

Angular momentum in the inner region (< 1 pc) is 
likely to be more complicated due to radiation pressure 
that would prevent collapse. The model of self-gravitating 
magnetized collapse in §^ was illustrative of the timescale 
involved but a detailed model for the inner region that 
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takes into account radiation pressure, radiative viscosity 
and evolution of the field is required. If dynamo action 
takes place within the compact region (< 100 pc), the ex- 
ponentiation timescale in the linear regime will be of or- 
der l/a;(ro) ~ 10^ yrs. The residual large scale seed field 
(~ lO^^G) from supernovae can be amplified to dynam- 
ical values (mG) in 10^ years and the e-folding timescale 
for magnetic braking is r^/vA ~ {Mir^) /rfj)^/'^ / B = 10^ 
yrs corresponding to this field strength. 



8. Conclusions 

In this work, a semi-analytic model of quasar formation 
was attempted which captures the essential details of star 
formation with supernovae feedback and angular momen- 
tum transport via magnetic fields and self-gravity. Based 
on the results in this work (as summarized in and 
for the relevant range in parameter space, it is seen that a 
IO^^Mq cloud at z of 8 can collapse to form a gaseous disk 
in the dark halo by z of 4.8 and about 10^ — 10* Af© of gas 
(which is 0.01-0.1 of the total disk mass) accretes into a 
'"0^ 100 pc region in a fast timescale of 10® yrs. The mag- 
netic stress from supernovae can be significant in the disk 
and is preferred over other possibilities as observational 
evidence exists for star formation, and large scale fields at 
z ~ 2. The collapse solutions do not hinge on the source 
of magnetic fields- SN or small scale dynamo. The small 
scale dynamo needs a few rotation time scales to build up 
and that could be significant- 10* yrs. The arguments for 



SN fields are given in section §4.1; the strong fields sug- 
gested can cause fast accretion rates instantly after the 
disk forms. The field strength estimation relies only on 
flux freezing arguments. Further, the alternative collapse 
solution using gravitational instabilty is also given. 

The solution to the non-linear diffusion equation pre- 
sented in §^ can be applied in other contexts, eg. proto- 
stellar disks. As a part of future work it is planned to 
investigate a detailed model of the inner region using a 
simulational approach, the impact of the restriction on 
the mass range of black hole hosts {M — z) by supernovae 
feedback on the details of the quasar luminosity function, 
dwarf galaxy formation, and the strength of the large scale 
seed field. 
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Appendix A: Core condensation in halos 

Here we calculate the change in gravitational potential 
energy from a truncated isothermal halo to a truncated 
isothermal halo of stars and dark matter containing an 
exponential gas disk. The final halo mass, = (1 — 
fg)M. is taken to be due to stars and dark matter, and 
remaining mass is in the form of a gaseous disk with Md = 



Mfb{l — /*) = fg-M where fg is the gas fraction. The 
initial potential energy is given by 



W,, = — 



2r„ 



(A.l) 



Ph^d) dV, (A.2) 



The final potential energy is given by 

2ry rd 2 J 

where the first two term represents the self potential en- 
ergies due a truncated isothermal halo of size r^, and an 
exponential disk of size rd, respectively, and the third term 
is due to the interaction between them. The first part 
of the last term with (the disk surface density S(r) = 
exp{-r/rd) Md/{2'Krj) and = GM?,/r«)can be re- 
duced to -{l/2){GMdMh/ry){l + exp{-l/ Xy)/ Xy) where 
the term involving A^, is negligible in the range of inter- 
est for which 0.1 > A„ > 0.05. The second part of the 
interaction term can be evaluated with the disk potential 
expressed as a Bessel series in the region r < where 
Ph = M/i/(47rr„r^) contributes. This term works out to 
be —{GMdMh/ry)0.3/Xv to a good approximation, where 
the c olla pse factor, r^/rd — 1/A^, is in the range 10-20 
(see |2.l| ). Combining the two we obtain 



^ J {Pd<^h + Ph^d) dV 



GMdMh 
2 r„ 



Finally, we obtain the change in gravitational potential 
energy 

GM^ 



/^W{MJg,Xv,Ty) 



2 r„ 
0.3 



(1 - fgY - 1 + _ + /,(1 - 



1 + 



GM^ 

2 r„ 



fa 



GM^ 



0.6 
A 

(U 1 

A„ 2 



/.h--o (A.3) 



which is bounded by 0.5 GM.^ /vy. 



Appendix B: Time-dependent evolution of an 
alpha disk 

Here we first consider a simple polytrope and then discuss 
the case of a disk coupled with an local energy dissipation 
condition. By assuming that the gas is a polytrope of index 

7, 

P^K {Y./2H)'^, (B.l) 
and using the vertical equilibrium, eqn (|3^, we obtain 
H^; = -2assPH 

= -a,,i^2/(7+l)s(37+l)/(7+l)^2(7-l)/(7+l). (B.2) 

In some applications it is appropriate to consider a disk 
with an energy dissipation condition as has been done by 
Cannizzo, Lee & Goodman 1990 (CLG), for the case of 
disk accretion of a tidally distrupted star onto a massive 
black hole. Using the solution presented in §||, we can gen- 
eralize to a case of a rotational law, a; oc r", and viscous 
stress, Hr0 cx S^'r'^. The viscous stress is given by 

H,^ = -2a..PH = -^IfMs. (B.3) 
pmp 
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where T refers to the central temperature and P is the gas 
pressure. Further, in the general prescription, the radiative 
flux is matched to the viscous dissipation by 

4 AaT'^ 

a...HP\a\ = - (B.4) 

where a is the Stefan-Boltzmann constant and T) is 
the opacity. This reduces to 



^3 ^ _3_ / assh 

16 \ /ITOp 



(B.5) 



Now, putting this back in eqn (B^), we obtain 



o \ 1/3 / , X 4/3 



16 



1/3 



|i/3 



(B.6) 



This implies b — 5/3, c — a/3 for a constant k. For the case 
of a Kepler potential we obtain the scaling (and opacity 
due to electron scattering which is a constant) used in 
CLG (a = -3/2, c = -1/2,6 = 5/3) and their solution 
is given by eqn (^) for these indices. The case of a flat 
rotation law and opacity due to electron scattering yields 
(a = —l,b = 5/3, c = —1/3). This solution can be easily 
extended to more general opacities of the form k cx p^T'', 
such as the Kramer's, and the appropriate solutions can 
be easily found. 
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